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In this paper we discuss the criticality of a quantum Ising spin chain with competing random 
ferromagnetic and antiferromagnetic couplings. Quantum fluctuations are introduced via random 
local transverse fields. First we consider the chain with couplings between first and second neighbors 
only and then generalize the study to a quantum analog of the Viana-Bray model, defined on a 
small world random lattice. We use the Dasgupta-Ma decimation technique, both analytically and 
numerically, and focus on the scaling of the lattice topology, whose determination is necessary to 
define any infinite disorder transition beyond the chain. In the first case, at the transition the 
model renormalizes towards the chain, with the infinite disorder fixed point described by Fisher. 
This corresponds to the irrelevance of the competition induced by the second neighbors couplings. As 
opposed to this case, this infinite disorder transition is found to be unstable towards the introduction 
of an arbitrary small density of long range couplings in the small world models. 
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Quantum fluctuations play a crucial role in the spin 
glass phases of the Sr-doped cuprate La2Cu04 llL or the 
dipolar glass LiHo^Yi_ a; F4 in a transverse field [2j. Ran- 
domly coupled quantum two level systems also appear in 
the understanding of the dielectric response of low tem- 
perature amorphous solids Q, and as the main low fre- 
quency source of decoherence of solid state quantum bits 
In all cases, the quantum fluctuations compete with 
the random couplings between the spin, and tend to dis- 
order the corresponding random ordered phases. 

One of the simplest random quantum model to study 
this competition is probably the random Ising spin model 
in a transverse magnetic field. 
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where af, af are the usual Pauli matrices, and the trans- 
verse fields hi are responsible for the quantum tunneling 
fluctuations between the up and down states of the Ising 
spins. In a pioneering work, Fisher has given asymptoti- 
cally exact results for this random quantum Ising model 
with first neighbors random ferromagnetic bonds in one 
dimension 0. By using a decimation technique devel- 
oped by Dasgupta and Ma 0, he described the infinite 
disorder quantum phase transition of this model. Some 
of the main features of this peculiar transition were a 
diverging dynamical exponent z, and very strong inho- 
mogencities manifesting through drastically different be- 
havior between average and typical correlation functions. 

Natural extensions of these results to higher dimen- 
sions have proved to be difficult. In particular, an analyt- 
ical implementation of the Ma-Dasgupta decimation be- 
yond the simple chain is extremely cumbersome. The rea- 
son is that any initial lattice except the chain is quickly 
randomized by the decimation. Thus one has to resort 



to a numerical implementation of this decimation 0, • 
For two-dimensional regular lattices, the results for the 
random ferromagnetic Ising model arc consistent with the 
survival of an infinite disorder quantum phase transition, 
albeit with exponents different from the one-dimensional 
case Q. On the other hand, the quantum Ising spin 
glass, corresponding to the model with both fer- 
romagnetic (positive) and antiferromagnetic (negative) 
couplings, was also studied in two and three dimensions 
via Monte-Carlo simulations 0. The numerical works 
found no sign of an infinite disorder quantum critical 
point. Since the results for the random ferromagnetic 
model are expected to extends to the quantum Ising spin 
glass, this discrepancy certainly deserves further work. 

In this perspective, we investigate in this paper the 
stability of infinite disorder fixed point of the quantum 
Ising spin glass chain with respect to competing further 
neighbors couplings in two extreme cases. In a first step, 
we focus in the case where second neighbors couplings 
are present in model (JTJ besides the first neighbors cou- 
plings. The couplings are taken as either ferromagnetic 
or antiferromagnetic. By combining analytical (for small 
second neighbors couplings) and numerical decimation 
techniques we investigate the relevance of the presence 
of higher range couplings and their interplay with ran- 
dom signs in the couplings. We pay a special attention 
to the topology of the renormalizcd lattice, which appears 
crucial in the precise characterization of infinite disorder 
transitions. 

A natural complement to this first case consists in con- 
sidering this quantum Ising spin glass on a random net- 
work, obtained by adding a finite density of long range 
couplings between the chain's sites. Indeed, our model 
can be considered as a quantum analog of the classical 
spin models of Viana and Bray [13, although we keep 
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a local regular topology besides the random long range 
couplings in our small world lattice^]]]. I n the case of 
classical spin glasses, these random lattice models are a 
natural extrapolation between the short-range model and 
its mean-field version. They undergo a finite temperature 
transition of the mean-field type, albeit with peculiarities 
induced by the finite connectivity Wether the infi- 
nite disorder physics survives to this tendency towards 
mean-field like physics is the natural question we will 
consider. 

In the first part of this letter, we consider the model (JTJ 
on a chain, with only first and second neighbors couplings 
(Zig-Zag ladder) |lj| . Both the first and second neighbors 

couplings Jii + i,Jii + 2 can be- antiferromagnetic (< 0) 
with probability p, and ferromagnetic with probability 
1 — p. The I-T^IliI are uniformly distributed between 

(21 (2) 

and 1, the |Jj between and JmLx, and the trans- 
verse fields hi between and h max . Note that via an 
appropriate unitary transformation we can map this sys- 
tem onto one where only the second neighbors couplings 
can have both signs, but at the cost of a modification of 
the magnetic properties of the system. Hence for clar- 
ity, we prefer to consider only the more natural choice 
defined above. 

We will analyze the low temperatures behavior of this 
system by means of the Dasgupta-Ma decimation tech- 
nique which was exploited by Fisher in the case, 
among others, of the random ferromagnetic quantum 
Ising chain. Its extension to the present case of mixed 
coupling (anti- ferromagnetic and ferromagnetic) contains 
one supplementary rule as detailed below. The running 
energy scale f2 is defined as the maximum of the ampli- 
tudes of bonds \J[j\ and fields hi. At each decimation 
step, if this maximum corresponds to a field hi, the cor- 
responding spin is frozen in the x direction, generating 
new couplings Jjk = Jjk + (JijJik/ty between all pairs 
(J, k) previously connected with the spin i. On the other 
hand, if the maximum is a ferromagnetic coupling J^, the 
two spins i and j are paired to form a new cluster [ij] 
of magnetization /Lt^-j = fa + \Xj (where \n corresponds 
to the magnetization of cluster i), and coupling with site 
& J[ij]k = Jik + Jjk M- The new rule occurs when this 
maximum coupling is anti-ferromagnetic. In this case, 
if e.g the magnetization pn is larger than then the 
new cluster's magnetization reads /ir^i = fii — fj,j, and 
the interaction with a third spin k is Jryifc = Jik — Jjk- 
In both cases, the effective transverse field acting on the 
new cluster is huj] = hihj/Q. 

An analytical study of the scaling behavior of the 
model Q under the above decimation rules is difficult 
even on the Zig-Zag ladder we consider. As mentioned 
in the introduction, couplings Jy are quickly generated 
on many length scales \i — j\, and the initial lattice is 
quickly randomized (sec below). To fix the notation and 
clarify the procedure, it is useful to start by consider- 



ing the evolution under the RG of the first neighbors 
chain, extending the result of Rcf. @ to the presence 
of anti-ferromagnetic couplings. We introduce the con- 
venient logarithmic variables pi = ln(0//ij), Ci,i+i = 
ln(0 / 1 t7j,i+i | ) and scaling parameter T := ln(f2o/^) 
where fio is the initial value of Q,. Their "distributions" 
densities are defined as TZ(/3, T) for the fields, V (1+) (C, H 
for the ferromagnetic bonds, and "PH(f, T) for the 
anti-ferromagnetic bonds. Note that while lZ(f3,T) is 
normalized, for the bonds only the sum T) = 

£>(!+) (C,r) + P( x -)(e,r) has a norm one. As can be 
deduced by a gauge transformation of Q), 1Z(/3,T) and 
7>W(C,r) satisfy the same differential scaling equations 
than in the ferromagnetic case Q provided we use the 
maximum instead of the sum in the above decimation 
rules, which is valid for broad enough distributions. Fi- 
nally, the function V{(3,T) = "P (1+) (C,r) _p(i-)(^r) 
is found to satisfy the same scaling equation than "PW . 
The fixed point TZ = = V*{x,Y) = e~ x / T /T of the 
ferromagnetic chainQ is easily extended to the two fol- 
lowing case : the above ferromagnetic point now corre- 
sponds to the solution V = V*, or V {1+) = = V*, 
pi 1 -) — 0. As expected, it can be explicitely checked in 
the RG equations that this fixed point is unstable towards 
the proliferation of anti-ferromagnetic bonds. The new 
transition point corresponds to the solution D = 0, or 
r p( 1 +) = = 7 ? */2, corresponding to an equal den- 

sity of random positive and negative couplings. Hence, 
we will loosely call it the spin glass fixed point by analogy 
with the physics of the classical model in higher dimen- 
sions. The characteristics of the transtion from the the 
ferromagnetic to the disordered phase obtained by Fisher 
[j| translate to the present spin-glass fixed point into an 
average linear susceptibility (under the application of a 
small z field h) which diverges as x(T) - |lnT|*- 2 / r 
where cf> = (1 + V5)/2 . Similarly, we extract the scal- 
ing behaviour of the average non linear susceptibility 
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ensemble average, as Xni(T) 

Having clearly defined the notation and fixed points 
for the chain, we can now study pcrturbativcly their 
stability with respect to small second neighbors compet- 
ing interactions. To first order, such an analysis can be 
conducted by considering the presence of J 1 - 2 ^ negligible 
compared to the and checking whether this condi- 
tion is self-consistcntly preserved under the re-scaling. 
More precisely, we will assume that (i) a J- i+2 will 
never constitute the highest energy in the system and 

(2) 

therefore never be decimated (ii) in sums, the Jl i+2 arc 

negligible with respect to JiJ+i (Hi) creation of third 
neighbour couplings out of second neighbour couplings 
can be neglected. As above, we define the distribution 
(C, r) := -p( 2+ > (C, r) +-P< 2 -) (C, r) as the sum of "dis- 
tributions" of positive and negative next nearest neigh- 
bour couplings. With the above hypothesis, its scaling 
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behavior is found to be described by 

^K) = ^K)_ pW(c) ( 2K(0) + P < 1)(0) ) 
+ 2K(0) / dCi d( 2 P (1) (Ci)^ (2) (C 2 )^(C-Ci -C2) 

+ V w (0)S(C-A) (2) 

where A is an arbitrary large constant which stands for 
the negligible J in log. coordinates, and is taken to 00 at 
the end of calculations. The T depcndance of the distri- 
bution has been omitted for clarity. With the above hy- 
pothesis, the probability distributions for fields and near- 
est neighbour couplings still follow the equations for the 
chain. Hence, at the "Spin Glass" critical point, we can 
insert the scaling form 1Z = V^ — V* in (J2J). It is usefull 

to split 7>( 2 )(C,r) into a A independent part vf ) {z,Y) 
(2) 

and 7-^ (z,r). By denoting p(z,T) the Laplace trans- 
form in £ of V( 2 \^,T), we find that for z and T finite 
and fixed, p\{z,T) — > when A — > 00. Then we show 
that the norm of the two parts of the solution satisfy : 

Unfile = 1 - Unfile = lim^ Pi(*,r) = r /r, corre- 
sponding to a constant "decrease" of the couplings J^ 2 \ 
In this regime, the system "forgets" its initial conditions 
and flows to a general state gouverned by . Con- 
sistency of conditions (i) to (Hi) can also easily being 
checked from the properties of the Laplace transform. 
As a consequence, such small next nearest neighbor cou- 
plings correspond to an "irrelevant perturbation" at this 
infinite disorder fixed point. 

To go beyond this perturbative analysis, we have stud- 
ied the scaling behavior of the zig-zag ladder by imple- 
menting numerically the above renormalization rules. We 
start by choosing a random configuration of fields hi and 
couplings i+D -Tj i+2 according to the previous initial 
distributions probabilities. Then at each step, the en- 
ergy scale is lowered and the number N of spins is re- 
duced by 1 according to the decimations rules specified 
above. This process is continued up to the last remain- 
ing spin, and repeated for a number R = 10 3 config- 
urations. No assumption is made on the topology of 
the renormalized lattice, and we keep a priori all gen- 
erated couplings. However, for practical reasons it ap- 
pears necessary to restrict ourselves to energies larger 
than a lower cut-off fl m in- With this procedure, the dis- 
tributions V(C,T),TZ(P^T) are correctly sampled below 
F m ax = hi(f2o/n m ; n ) I?) . For most of our results, this 
cut-off fl m in was maintained to negligible values, with- 
out any noticeable incidence on the results. For fixed 

(2) 

Jmax-i the transition is reached by varying the maximum 
amplitude h max of the fields. We locate a putative infi- 
nite disorder phase transition by using the analogy with 
percolation At each decimation step i, correspond- 
ing to a system size No — i, we consider the number of 
realizations rih{i) where a field was decimated at step i, 




c : number of neighbours per spin (connectivity) 

FIG. 1: Scaling behavior of the distribution of connectivity 
P(c) for different number TV of remaining clusters (spins), for 
the Ji — J2 model. The initial size is iVo = 2 14 = 16384 spins, 
and the decimation was performed over 1000 samples. After 
a transient regime characterized by an algebraic distribution 
of connectivity, the distribution ultimately renormalizes to- 
wards a delta function c = 2 corresponding to the topology of 
the chain. The inset show the ratio of number of decimated 
fields over the number of decimated bonds as a function of 
the number of remaining clusters N (see text). The critical 
point is evaluated as h max = 2.20 ± 0.05. 



the similarly for the bonds nj{i). At the transition, the 
ratio rih(i)/nj(i) should become scale invariant, whereas 
it should diverge or decrease to zero respectively in the 
disordered or ordered random phases. Moreover, the scal- 
ing behavior of this ratio is an excellent way to check for 
possible finite size effects respective to the topology of the 
initial lattice. The inset of the figurc^shows this scaling 
of the decimation ratio for two values around the candi- 
date critical value of h max . Once such candidates for the 
transition are determined, we have studied the scaling 
behavior of the distributions functions P((,T),TZ((3,T), 
of the distribution of magnetization ^t(F), and number of 
active spins n(T) in the clusters. This allows to charac- 
terize the criticality of the infinite disorder fixed point. 
Moreover, to fully characterize an infinite disorder fixed 
point beyond the simple chain, one should also be able to 
determine the renormalized topology of the critical lat- 
tice, and the associated correlations with the couplings. 
In a first attempt to study the scaling of this topology, we 
have followed the distribution of the connectivity of the 
lattice as the decimation goes on. The results, depicted 
on fig. shows that while initially all sites have only 
4 neighbors, the distribution P(c) flows towards an in- 
termediate algebraic distribution at intermediates sizes. 
While highly connected sites appear, we find by varying 
our lower cut-off T max that rather strong correlations ex- 
ist between the bonds connecting these sites. And while 
the decimation is pursued, the distribution narrow back 
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towards a delta function peaked on c = 2, i.e. the lat- 
tice is ultimately rcnormalizcd towards a chain. We thus 
find that for the J\ — J2 model, the infinite disorder fixed 
point is always given by the fixed point of the chain (see 
above and [^| ) , in agreement with previous results on the 
similar ferromagnetic two-leg ladder Q 




Number of remaining spins 



ladder) leads to an asymptotic lattice equivalent to a sim- 
ple chain, proving the irrelevance of the second neighbor 
couplings perturbation at the infinite disorder fixed point 
of the chain. On the other hand, the results of our nu- 
merical renormalization approach show that the inclusion 
of an arbitrary density of long range couplings in the 
chain modifies the scaling behavior of the lattice's topol- 
ogy, and thus the associated critical behavior. These re- 
sults stress the importance of determining the renormal- 
ized topological properties at any possible infinite disor- 
der transition beyond the one-dimensional examples. In 
particular, the intermediate regime we have identified in 
our study of the Zig-Zag ladder opens the possibility of 
new infinite disorder scenarii for models with correlated 
long-range couplings. A natural extension of the present 
work would certainly focus on random algebraic interac- 
tions and the effect of the dimension, possibly relevant 
to the understanding of the dipolar glass LiHo 2 Yi_ l: F4 
in a transverse field |2| • 

D. Carpentier and P. Pujol would like to acknowledge 
F. Leonforte for collaboration in a preliminary investiga- 
tion related to the present work. 



FIG. 2: Decimation ratio for the disordered quantum Ising 
model on a small world lattice. The initial size is No = 1024, 
q = 0.1 and 5000 samples were used. The scaling behaviour of 
the distribution of connectivity for h max = 1.5, in the inset, 
shows a broadening of this distribution up to some clear finite 
size topological effects inherent to small world models. 

The previous results motivated the study of the op- 
posite limit of long-range couplings competing with the 
initial couplings of the chain. Thus we naturally consider 
the hamiltonian JTJ on a small world lattice[ll|, where 
beyond the previous nearest neighbors couplings 
we add random infinite range couplings J/^ between any 
two non-neighbor sites i and j, with density p/N. In 
this paper, the existing couplings J/^ and J-i +1 are dis- 
tributed with the same uniform distribution between ad 
1. With these conventions, the average initial connectiv- 
ity of this lattice is 2 + q. Results of the same numerical 
decimation procedure as above indicate a phase transi- 
tion different from the previous one (Zig-Zag ladder). In 
particular, contrarily to the previous case, the distribu- 
tion of connectivity of the rcnormalizcd lattice broadens 
without limit up to some finite size effects. This cross- 
over happens when the numerical upper bound of the 
renormalized distribution P(c) becomes of the order of 
the system size. Once this happens, highly connected 
sites proliferate, leading to a mean-field like behavior. 

In this paper we have shown how the presence of ran- 
dom signs and further neighbor couplings affect the crit- 
ical behavior of the random quantum Ising chain. We 
have particularly focused on the topological properties 
of the renormalized lattice, and we have explicitly shown 
how the presence of second neighbors couplings (Zig-Zag 
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